Strategies for the preparation of large cluster states using non-deterministic gates 
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The cluster state model for quantum computation has paved the way for schemes that allow 
scalable quantum computing, even when using non-deterministic quantum gates. Here the initial step 
is to prepare a large entangled state using non-deterministic gates. A key question in this context 
is the relative efficiencies of different 'strategies', i.e. in what order should the non- deterministic 
gates be applied, in order to maximize the size of the resulting cluster states? In this paper we 
consider this issue in the context of 'large' cluster states. Specifically, we assume an unlimited 
resource of qubits and ask what the steady state rate at which 'large' clusters are prepared from this 
resource is, given an entangling gate with particular characteristics. We measure this rate in terms 
of the number of entangling gate operations that are applied. Our approach works for a variety of 
different entangling gate types, with arbitrary failure probability. Our results indicate that strategies 
whereby one preferentially bonds together identical qubits are considerably more efficient than those 
in which one does not. Additionally, compared to earlier analytic results, our numerical study offers 
substantially improved resource scaling. 
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states (cluster states in the case of Ref. [2.] and 'linked 
states' in Ref. [3]) using non-deterministic linear opti- 
cal gates, and calculate the efficiency of this process for 
particular values of the success probability of the elemen- 
tary entangling operations, Pgate, with Pgate > 1/2. Sub- 
sequently, in Ref. 0], a simplified linear optical scheme 
for generating cluster states was described, and it was 
also shown how to scalably construct clusters in the case 
where Pgate = 1/2. In Ref. ^8|] a comparable scheme is 
described within the standard circuit model. 

In Ref. Q , and subsequently Ref. , these ideas where 
extended to the case of arbitrary values of Pgate- In partic- 
ular, it was shown that, even for arbitrarily small values 
of 

Pgate? linear clusters can be grown with a cost that 
grows linearly with the size of the required cluster, al- 
beit with a linear coefficient that becomes rather large 
for very small values of Pgate • 

One shortcoming of the calculations presented in 
Refs. 0, is that, in order to obtain analytic approx- 
imations for this overhead cost, certain assumptions are 
made that necessarily reduce the efficiency of the pro- 
tocol. In particular, it is assumed that short clusters 
(below some critical length that depends on Pgate) are 
constructed without recycling - upon gate failure, the 
remaining sub-clusters are discarded and construction 
starts from scratch. In this paper, we avoid such prof- 
ligate measures by directly computing the cost of a more 
frugal scheme via numerical simulations. 

In general, the questions of optimal strategies and 
physical resource requirements are difficult problems. 
Previous work by Kieling et al. [TO] considered this 
issue. Here the authors specifically focussed on the lin- 

ear optics implementation of Ref. [7] , where gate success 

probability is assumed to be Pgate = 1/2. They assumed 
1rohde@p hysics . uq.edu . au URL: ht tp : //www . physics . uq . edu . au/peop£fe^ :egKa^p^ce of N Bell pairs, and asked, for a given strat- 



I. INTRODUCTION 

The development of cluster state (or one-way) quan- 
tum computation [l| fias presented us with a completely 
alternate perspective to the standard circuit model for 
quantum computing (QC). This model has the benefit 
that it transfers all entangling quantum gates to an of- 
fline state preparation stage. This has attracted much 
interest in the quantum computing community because 
in some instances it has highly favourable properties. One 
of the key benefits offered by this model is its applica- 
bility to architectures where entangling gates are non- 
deterministic. Here the cluster state model can allow sig- 
nificant physical resource savings when compared to the 
circuit model, most notably in optical implementations 
[I, H, 0, m, 0, . The central idea here is that a large clus- 
ter state, a highly entangled state of many qubits, can be 
'grown' using a 'divide-and-conquer' approach, whereby 
smaller clusters are iteratively bonded together, a process 
which can lead to very efficient physical resource scaling 
characteristics. While several approaches for efficiently 
growing cluster states using non-deterministic gates have 
been prescribed, a key open question is "what are the 
minimal physical resource requirements for preparing a 
certain cluster using entangling gates with a given success 
probability?" and "what strategy (or algorithm) must be 
employed to achieve optimal resource scaling?" . 

A number of other authors have considered the prob- 
lem of efficiently constructing cluster states using non- 
deterministic two-qubit entangling operations. Refs. pl,[3| 
present procedures for generating multi-qubit entangled 
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egy of bonding them, what the expected length of the 
longest chain, Q{N)^ will be. They then considered the 
behavior of Q{N) against N for various strategies and 
asked what the upper bound on this scaling relationship 
was. There are two strategies that were given special at- 
tention in their analysis that we will repeat here. The 
first is Greed, whereby we always preferentially bond 
the largest existing cluster states together. Intuitively one 
might expect such a strategy to perform extremely well 
- we rapidly build up clusters of higher length. The sec- 
ond is Modesty, where we do exactly the opposite and 
always preferentially bond the smallest available clusters 
together. Surprisingly, the authors showed that Greed is 
in fact not only a sub-optimal strategy, but MODESTY is 
substantially better and close to optimal. More recently, 
other strategies for growing cluster states based on per- 
colation theory have been examined pj| . 

In this paper we consider the question of cluster state 
preparation strategies from a slightly different perspec- 
tive. Specifically, if we are ever to implement large scale 
cluster state quantum computing we will effectively be 
operating in a regime where we must have an effectively 
infinite resource of qubits. In such a situation we are in- 
terested in the rate at which we can prepare 'large' clus- 
ters from this resource, i.e. how many elementary entan- 
gling operations must we perform per large cluster we 
prepare? Thus, in our analysis we assume an infinite re- 
source of qubits and calculate the rate at which large 
clusters are prepared from this resource, given a particu- 
lar entangling gate. Our approach is to employ a numeri- 
cal Monte-Carlo simulation of the preparation procedure. 
This technique works for a variety of different entangling 
gate types with arbitrary gate success probability. 

In this paper we focus on the preparation of linear 
clusters, which are basic primitive in the construction 
of larger clusters. There are many strategies for prepar- 
ing arbitrary large cluster states. For example, in the 
scheme of Nielsen [2], a resource of star-shaped 'micro- 
cluster' is used to prepare large cluster states using non- 
deterministic gates. These star-shaped clusters can in 
turn be prepared from linear clusters using local oper- 
ations and single qubit measurements. Similarly, in the 
scheme of Duan & Raussendorf [6| a resource of linear 
clusters are used to prepare '-|- '-shaped clusters, which 
facilitate the construction of lattice clusters using_^non- 
deterministic gates. The scheme of Kieling et al. [9|, [lo| 
describes a 'weaving' technique, where a lattice cluster is 
prepared directly from a number of linear clusters. A dif- 
ferent approach to preparing lattice clusters from linear 
clusters is described in Ref. [4, 5, 7, 12]. In these cases the 
additional cost of preparing 2D clusters from ID chains 
is well understood. Therefore, it is appropriate to focus 
our attention on the preparation of linear clusters as a 
primitive for use in such higher- level protocols. 

This paper is structured as follows. In Section [III we 
present background material on cluster states, strategies 
and gate types. In Section IIIII we describe our analysis 
techniques, and in Section [iVl we discuss our simulation 



results. We conclude in Section [Vl 



II. BACKGROUND 

A. Cluster states 

The standard circuit model for quantum computation 
is very analogous to our usual understanding of classical 
circuits - we begin with an input state, apply a series 
of gates to manipulate the state, and then measure the 
output to determine the 'answer'. The cluster state model 
is completely different and has no analogy in classical 
circuit theory. Here we prepare a highly entangled state, 
known as a cluster (or graph) state, and then perform 
some sequence of single qubit measurements. The choice 
of measurements - both the measurement basis and the 
order in which they are performed - then implements the 
algorithm. 

It has been shown that lattice clusters are universal for 
quantum computations up to a given size. Thus lattice 
clusters act as a resource for universal quantum compu- 
tation. For this reason theoretical studies often focus on 
lattice clusters. 

A cluster state can be expressed as a graph, hence the 
alternate term 'graph state'. In this representation ver- 
tices represent single qubits that are initially prepared 
in the |+) = (|0) + |l))/>/2 state. Graph edges repre- 
sent controlled-phase (CZ) gates performed between the 
neighbouring vertices. This bonding procedure is shown 
in Fig. [H 
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FIG. 1: Creating links in a cluster state using a CZ gate (top). 
Removing a qubit from a cluster state with a Z measurement 
(bottom). 

An additional property of cluster states that we will 
make use of in this paper is qubit removal. When a qubit 
is measured in the Z eigenbasis this removes the qubit 
and any edges connecting to it, while leaving the remain- 
der of the cluster unchanged. This property is useful to 
recover from certain types of gate failure. For example, 
consider a linear cluster to which we attempt to bond an- 
other cluster at the end. If the bonding fails with no well- 
defined failure outcome (i.e. the qubit is effectively depo- 
larized) then performing Z measurements on the qubits 
adjacent to the affected ones recovers the remainder of 
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Strategy 


Choice of clusters to bond 


Greed 
Modesty 
Random 
Paired Greed 
Paired Modesty 
Paired Random 


Two largest available clusters 
Two smallest available clusters 
Randomly choose two clusters 
Largest available identical clusters 
Smallest available identical clusters 
Random identical clusters 



TABLE L Summary of different bonding strategies. For the 
Modesty and Paired Modesty strategies single qubits are 
excluded from the decision making process, unless no larger 
clusters are available. This is necessary since we assume an in- 
finite resource of single qubits. This is distinct from the Mod- 
esty strategy employed in Ref . [1, |lp] . The Paired variations 
of strategies always only bond two identical clusters, and bins 
containing only a single chain will undergo bonding. 



the cluster by isolating the malignant qubits from the 
rest of the graph. This is shown in Fig. [TJ 

We have only recapped the basic properties of clus- 
ter states requisite for our analysis techniques. Excel- 
lent reviews of the cluster state model can be found in 

Kefs. 



B. strategies 

Beginning with a resource of qubits, we begin bonding 
cluster together to construct progressively larger cluster 
states. At any point in time we will have a 'pool' of clus- 
ters of varying lengths, to which we may apply bonding 
operations. A 'strategy' is an algorithm that inspects the 
current pool of available clusters and decides which two 
clusters to attempt to bond together. Previously we re- 
viewed two strategies - Greed and Modesty. In our 
simulations we will consider several additional strategies, 
which we summarize in Table |I] (Note that the Mod- 
esty strategy applied in this work differs from that in 
Ref. as will be discussed in Section [IV|) . 

Typically there will be much room for parallelizabil- 
ity in bonding strategies, limited only by the number of 
gates available, the number of clusters in the pool, and 
potentially other practical constraints. While this paral- 
lelisation will not reduce the overall number of bonding 
operations that must be performed, it will, in general, 
reduce the average amount of time qubits spend in quan- 
tum memory. Clearly this will generally be highly favor- 
able from a decoherence perspective. 

In all our simulations do not employ any paralleliza- 
tion. That is, exactly one gate operation is applied per 
time step. There are two reasons for this. First, it elimi- 
nates a parameter and makes the results much easier to 
digest. Second, the resource efficiency of a strategy with- 
out parallelization will represent an upper bound on the 
achievable resource efficiency with parallelization, and 
therefore is an extremely useful parameter to know. This 
is because any strategy involving parallelization can be 



simulated without parallelization, but the converse is not 
true in general. Therefore, parallelizability represents an 
added constraint on the class of strategies that can effec- 
tively be implemented. 



C. Gate types 

In addition to various strategies for bonding clusters 
together there are a variety of gates that may be em- 
ployed to perform the necessary bonding operations to 
prepare cluster states. Specifically, there are two features 
that differentiate bonding gates. The first is their success 
probability, Pgate- In our simulations we treat this as a 
parameter. The second is their resource usage. That is, 
how qubits are wasted upon success and failure respec- 
tively. Formally, a bonding operation will act on a set 
of two linear clusters of lengths li and which, upon 
success or failure will be mapped to a set of new linear 
clusters. In Table [2] we summarize the most well known 
gate types and their resource usage. The action of these 
gates is shown graphically in Fig. [31 

In our simulations we will primarily focus on three 
gates in particular: two variations of the CZ gate - the 
archetypical bonding gate; and, the so-called EO (Entan- 
gling Operation) gate introduced in Ref. [4]. These gates 
are particularly relevant to current architectures for per- 
forming quantum computation using non-deterministic 
gates. 

The regular CZ gate, upon success, acts in the usual 
way, implementing a non-destructive controlled-phase 
operation. Upon failure its action is undefined. Thus 
to recover the remaining cluster one must measure the 
qubits neighboring the ones acted upon by the CZ gate 
in the Z-basis. This disconnects those qubits from the 
cluster resulting in two clusters, each two qubits smaller 
in size [15]. 

The KLM CZ gate differs in that upon failure it 
performs Z-measurements on the affected qubits. Thus, 
these qubits are automatically removed from the clusters 
and no additional measurements are needed to recover. 
Therefore, the KLM CZ gate is more resource savvy and, 
upon failure, results in two clusters each with only one 
qubit removed. 

The EO gate is employed in a hybrid scheme, involv- 
ing both photons and 'matter qubit' systems. These may 
be trapped ions or atomic systems, impurities in semi- 
conductors or insulators (e.g the NV-diamond system) 
or quantum dots. In this approach, the matter systems 
contain additional optical transitions, such that photons 
can be emitted whose state is conditional on the logical 
qubit state. In Ref. [4] the photons are encoded in sep- 
arate time windows, although other encodings such as 
frequency have subsequently been proposed [15]. 

Two matter qubit systems can be entangled by first 
inducing them to emit photons, passing these photons 
through a simple passive linear optical device (typically 
a beam splitter), and detecting the output modes with 



4 



destructive photo-detectors. In Ref. two such rounds 
are implemented, one for each time window. Success of 
the scheme is 'heralded' by a particular sequence of de- 
tector clicks, and results in a highly entangled state of 
the two qubits. 

In fact, it can be shown that the scheme amounts to a 
partial parity measurement. That is, in the case of suc- 
cess, a projection operator 11+ = |01)(01| + |10)(10| (up 
to local unitary corrections) is applied to the state of the 
two qubits. In the case of failure, each of the two qubits 
may be affected by a Pauli Z (phase flip) error. Success 
occurs with probability Pgate = where 77 is the com- 

bined photon collection and detection efficiency for each 
of the modes output from the beam splitter. 

This partial parity measurement can be used to gener- 
ate linear cluster states as follows. Pairs of unent angled 
qubits are first prepared in the product state |+)|+). Ap- 
plying the EO, followed by single qubit local unitary cor- 
rections [4] leads to a two qubit cluster state in the case 
of success. In the case of failure, the matter qubits are in 
an unknown state and should be re-initialised. Given an 
/i-qubit linear cluster state, a single qubit can be added 
to one end of the chain by preparing the qubit in the 
state 1+) and applying the EO between this qubit and 
the qubit at the end of the chain. On success, this results 
(up to local unitary corrections) in a chain of length /i + l. 
On failure, the qubit at the end of the linear cluster may 
have undergone a Z-error, and so this qubit should be 
removed from the cluster by measuring it in the Z-basis. 
The result is a chain of length — \ and two unent angled 
qubits which may be re-prepared. Finally, in this paper 
we will often consider the case when two chains of arbi- 
trary lengths, l\ and ^2, are fused at their endpoints. In 
this case, on success, the result (up to local unitary cor- 
rections) is a chain with a 'backbone' of length /i + ^2 — 1 
with a further dangling node (or 'cherry') which is an ex- 
tra qubit connected to the backbone at the point where 
the two chains were connected [12,]. For the purposes of 
this work, we assume for simplicity of our calculations 
that this cherry is removed via a Z- measurement. How- 
ever, in general the cherry may be utilized to increase 
the efficiency of growing linear clusters or higher dimen- 
sional graph states [l?]. Thus the estimates of the cost 
presented in this paper may be reduced even further. 

On failure, again the end qubits of each chain must 
be measured in the Z-basis, and the result is two chains 
of length l\ — \ and ^2 — 1, and two unentangled qubits. 
These results are summarized in Table [2j Note that, apart 
from the initial step of making two qubit clusters, the 
resource usage rules for the EO and the type-I fusion 
are identical. 



III. ANALYSIS 

We now turn our attention to our analysis technique. 
Our analysis proceeds as follows. We begin with a pool of 
linear clusters of various lengths. We describe the state 
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FIG. 2: Resource usage of different cluster state bonding op- 
erations, acting on the initial cluster states {h^h}- 'GZ' is 
a generic non- destructive cont rolled-sign gate. No qubits are 
lost upon success. However, the failure mode of the gate is 
undefined. Thus, recovering the clusters upon failure requires 
measuring the neighbouring qubits in the Z eigenbasis 0. 
The 'KLM GZ' gate is a Knill-Laflamme-Milburn [16] type 
controlled- phase gate which performs Z measurements on 
both qubits upon failure. This has the effect of removing the 
qubits from the cluster and no further recovery measurements 
are necessary. The 'type-I fusion' [7] and 'EO' [4] gates im- 
plement a partially destructive fusion. Upon success one of 
the two qubits is consumed, and upon failure both qubits are 
effectively measured in the Z eigenbasis. The 'type-II fusion' 
gate differs only in that upon success both qubits are con- 
sumed. 



of the pool by a population vector n, whose elements 
correspond to the number of clusters in the pool of the 
respective length. Beginning with an initial population 
vector no, we repeatedly apply a bonding strategy, which 
inspects the population vector and decides which two 
chains to attempt to bond together using the given gate 
operation. Let fs represent one iteration of the strategy. 
Then the evolution of the pool of clusters takes the form, 

fit+i = fs{nt)- (1) 

As discussed previously, we always assume strategies ap- 
ply exactly one gate operation per time unit. Thus, t gives 
the total number of gate operations that are applied. 

In our simulations we assume the population vector is 
initialized such that all elements are initially zero, except 
the first element (i.e. / = 1, the single qubit bin), which 
is set to infinity. As we begin applying a given strategy 
to the initial population vector, higher length bins be- 
gin to populate, and eventually longer chains will begin 
'spilling over' off the end of the vector. We regard any of 
these spillover chains as being 'complete'. Note that the 
definition of 'complete' depends on how many bins are 
used in the simulation, which we label L. We tally the 
total number of qubits that spill out of the population 
vector, which we label Nq. 

After a sufficient number of applications of the bonding 
strategy the system reaches a quasi steady-state, whereby 
the average population of each bin in the population vec- 
tor remains constant, averaged over a sufficiently large 
time period [17]. In our simulations we are interested in 
the quasi steady-state dynamics of this evolution. In par- 
ticular we are interested in the rate at which large clusters 
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FIG. 3: Bonding of linear clusters using the CZ, KLM CZ, 
type-I fusion, type-II fusion, and EO gates. 

are prepared per gate operation. Thus we define, 

r=limM. (2) 

t^OO t 

Note that 1/r corresponds to the average number of at- 
tempted entangling gate operation per qubit added to the 
completed chains. At first our treatment of large clusters 
may appear somewhat dubious, in that we are treating 
all large clusters, irrespective or their individual lengths, 
as being equal and simply putting them in one basket. 
The key observation here is that the resource overhead 
in bonding together multiple large clusters is negligible 
compared to the size of the cluster. Specifically, with a 
gate success probability of Pgate one expects to perform 
1/Pgate bonding operations on average until two clusters 
are successfully joined. If the lengths of the individual 
chains are much larger than this, /i,2 ^ l/Pgate^ this over- 
head is negligible and bonding can essentially be treated 



as 'free'. Thus, in the large cluster regime, it is justified to 
assume that a pool of large clusters is equivalent to a sin- 
gle large cluster of size equal to the cumulative size of the 
pool. Furthermore, for L 1/pgate, ^ becomes indepen- 
dent of L. Thus, for simulation purposes we can simply 
choose some value of L, sufficiently large that changes in 
L leave r unchanged. 

It is important to note that this technique is only 
suited to the analysis of linear cluster state preparation. 
The reason for this is that the population vector, n, con- 
tains a single entry for each possible graph configuration. 
When restricted to linear graphs the dimensionality of n 
is just the maximum number of qubit s. However, in the 
2D case the dimensionality of n grows exponentially with 
the number of qubits, ruling out even moderately large 
simulations. 



IV. RESULTS 

We simulate the effect of bonding strategies and gate 
parameters using a Monte-Carlo approach. We use a pop- 
ulation vector with L = 50 bins, and simulate over 50,000 
time steps. At each time step we tally Nq, the number 
of spillover qubits. Finally, at the end of the simulation 
we calculate r, the rate at which large clusters were pre- 
pared. We consider all six strategies listed in Table [D 
There is one caveat associated with the Modesty and 
Paired Modesty strategies. Specifically, we exclude the 
single qubit bin from the calculation, unless no other bins 
are populated. The reason for this variation, relative to 
the Modesty strategy discussed in Ref. [9, 10], is that we 
are assuming an infinite resource of single qubits. Thus, 
without this modification Modesty strategies would for- 
ever be stuck bonding single qubits together. 

In Fig. m we plot r against Pgate for all strategies. There 
are several interesting features. In agreement with [TO] 
the Modesty strategy performs better than Greed for 
all Pgate- The Random strategy performs worse again. 
Particularly interesting are the three Paired strate- 
gies, which all exhibit identical performance within the 
accuracy of our simulation. Perhaps surprisingly the 
Paired strategies perform significantly better than the 
un- Paired ones. In particular, for small values of Pgate, 
r drops to zero (within the accuracy our simulations) for 
all the uu-Paired strategies, which is not the case for 
the Paired strategies. 

We have compared the performance of a variety of 
strategies for a given gate, the CZ gate. Next we consider 
the performance of different gates for a fixed strategy. In 
Fig. [5] we plot the performance of the Paired Greed 
strategy for three different gate types: CZ, KLM CZ and 
EO. 

Consider the general behavior of these plots. In the 
limit Pgate = 1 the rate of large cluster preparation is 
r = 1. This is expected. It tells us that once the system 
attains its quasi-steady state, large clusters are spilling 
off the end of the population vector at an average rate 
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FIG. 4: Spillover rate against gate success probability for all 
six bonding strategies. The population vector has size L — 
50, and we apply a Monte- Carlo simulation over 50,000 time 
steps (5 X 10^ for lower plot). The lower plot is zoomed in 
on the region pgate < 0.2. All three Paired strategies exhibit 
identical performance, which has been labeled Paired. 



of one qubit per gate operation, which clearly must nec- 
essarily be the case. For Pgate < 1 we observe a lower r. 
This is because now the gate has some failure probabil- 
ity which results in larger clusters being converted into 
smaller cluster with non-zero probability. In other words, 
the dynamics of the population vector are no longer de- 
scribed by a one-way coupling from lower bins to higher 
bins, but now includes reverse couplings which slow down 
the overall flux of qubits towards higher length. 

The most striking feature of these plots is that all three 
Paired strategies exhibit identical performance, which 
is signiflcantly better than any of the other un-PAlRED 
strategies. Although we considered only a limited set of 
strategies, empirically this suggests that one should al- 
ways preferentially bond clusters of equal length, and be- 
yond this the performance is strategy independent. 



paired greed 




0.2 0.4 0.6 0.8 1 

Pgate 



FIG. 5: Spillover rate against gate success probability for three 
different gate types - CZ, KLM CZ and EO. All plots are using 
the Paired Greed strategy. The population vector has size 
L = 40, and we apply a Monte- Carlo simulation over 100,000 
time steps. Note that in practise the EO gate has a maximum 
success probability of p = 1/2. 



A. Comparison with analytic results 

We flnish by comparing the scaling relationships ob- 
tained through our numerical methods with two recent 
analytical studies of resource scaling. In the flrst in- 
stance we compare against the analytic results of Duan 
& Raussendorf [g^ for preparing linear clusters using non- 
deterministic CZ gates. In the second instance we com- 
pare against the analytic expressions obtained by Barrett 
& Kok for preparing linear clusters using the EO gate. 
See Appendix [A] for a discussion on the derivation of the 
analytic expressions. Note that in both plots the best 
numerical scaling relationships are signiflcantly more effl- 
cient than the analytic ones, owing to the lack of recycling 
in the analytic derivations leading to inefflcient resource 
usage. In the case of the EO gate we have included a third 
strategy, labeled EO Greed Paired, which is identical 
to the usual Greed Paired strategy except that if a Bell 
pair is available we always preferentially bond it with a 
single qubit to attempt to form a GHZ state. The reason 
for this exception is that with this type of gate bond- 
ing two Bell pairs together also forms a GHZ state upon 
success. Thus, doing so wastes the operation required to 
form one of the Bell states. At higher levels this discrep- 
ancy no longer occurs, so the usual pairing approach is 
favorable. 



V. CONCLUSION 

We have considered the question of preparing large 
cluster states using non-deterministic and potentially de- 
structive bonding operations, and the effect that bond- 
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FIG. 6: (top) Numerical scaling relationship for the CZ gate 
compared with the analytic scaling relationship obtained by 
Duan & Raussendorf. (bottom) Two numerical scaling rela- 
tionships for the EO gate compared with the analytic scaling 
relationship obtained by Barrett & Kok. The fist numerical 
simulation is for the normal Greed Paired strategy. The sec- 
ond, EO Greed Paired, is a modification of this where we 
always preferentially bond {h = 2,/2 = 1} pairs. 



ing strategies have on the resource requirements of this 
process. Empirically, our results indicate that Paired 
strategies, i.e. ones where we only bond clusters of equal 
length, perform much better than un- PAIRED ones, espe- 
cially in the regime of small Pgate • Additionally, all three 
Paired strategies considered exhibited identical perfor- 
mance. This suggests that the performance of large clus- 
ter state preparation is strategy independent provided 
one only bonds chains of equal size. This is in stark 
contrast to un-PAlRED strategies, where performance is 
highly strategy dependent. In this case it is observed that 
the Modesty strategy is superior to a Greed approach, 
consistent with the observations made in Ref. [§, [l^ . 

We also compared the Paired strategies with earlier 
analytic estimates for the cost of preparing linear clus- 
ters. We found substantial improvement in efficiency over 



these earlier estimates. We attribute this to the use of full 
recycling in our numerical analysis, whereas previous an- 
alytic results did not make full use of cluster recycling. 
This suggests that numerical techniques play an impor- 
tant role in the evaluation of resource scaling efficiency. 

The analysis techniques we described are very general 
and hold for arbitrary gate success probabilities and a 
variety of different entangling gate types, including de- 
structive ones. Therefore, our approach is suited to cal- 
culating physical resource requirements in a diverse array 
of different situations. However, our technique does have 
the disadvantage that it is only well suited to directly 
analyzing the preparation of linear clusters. 

Importantly, in our analysis we have not proven which 
strategies are optimal, we have simply made empirical ob- 
servations that some classes of strategies are better than 
others, and that some are equivalent. Thus there is much 
room for further analysis in this direction. The MAT- 
LAB source code used for our simulations is available at 



http: //www. physics .uq. edu. au/people/rohde/ 
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APPENDIX A: ANALYTIC ESTIMATES OF THE 
COST OF GROWING CLUSTER STATES WITH 
NON-DETERMINISTIC OPERATIONS 

An analytic estimate of the cost of growing cluster 
states using non-deterministic operations was given in 
Ref. 0. For the purposes of obtaining simple expressions 
for the cost of making linear clusters, the method used 
was intrinsically less efficient than the one used in our 
numerical simulations in this work. Here, we elaborate 
on the calculation presented in Ref. [4] , for the purposes 
of a comparison with our numerical results. 

The method used was to build short cluster chains 
(above some critical length, as described below) using an 
inefficient technique, and then join these chains together 
to form a long chain. Assume that we already have a lin- 
ear chain of length TV, and would like to add a short chain 
of length m to one end of this chain, using the EO with 
success probability Pgate- Using the rules given in Table 
[21 the expected length of the chain after this operation is 



L = PgM^ + m - 1) + (1 - Pgate)(A^ - 1) 

= N ^ mpgate - 1- 



(Al) 



8 



Thus, in order for the chain to grow on average, we re- 
quire L > N which imphes that the short chains must 
satisfy 



m > 



1 



Pgate 



(A2) 



Therefore the critical length depends directly on the gate 
success probability. 

To make these short chains of length m, we use a 'di- 
vide and conquer' technique. For ease of calculation, we 
assume that if any of the gates fail, any cluster chain frag- 
ments are simply discarded, and the process starts from 
scratch. For instance, to make a chain of length m = 3 
qubits, we attempt to bond 2-chains with single qubits. 
The total number of attempted gate operations to make 
a single 3-chain is then 



R2 

Pgate 
1 



1 



Pgate 
1 



Pgate Pgate 



(A3) 
(A4) 



Here, R2 = 1/pgate is the number of attempted opera- 
tions required to make a 2-chain. Similarly, we have for 
4-chains 



-R4 = 



1 



1 



Pgate 



Pgate 



Pgate 



(A5) 



To make longer chains using this inefficient procedure, 
we can recursively join these shorter chains together. For 
instance, to make a 5-chain, two 3-chains can be joined 
together. Joining two 5-chains together results in a 9- 
chain, joining two 9-chains results in a 17-chain, and so 
on. If all this is done without recycling, the number of 
attempted operations to make 5-, 9-, and 17-chains, re- 
spectively, are given by 



i?5 

Rg 
Rl7 



2 


2 


1 


3 

Pgate 

4 


^ 2 - 

Pgate 

4 


Pgate 

2 


Pgate ^ 

8 


^ 3 - 

Pgate 

8 


^ 2 

Pgate 

4 


Pgate 


Pgate ^ 


^ 3 

Pgate 



Pgate 

2 

Pgate 



Pgate 



(A6) 



One can easily recurse these relations further to deter- 
mine Rm for larger values of m (corresponding to smaller 
values of Pgate). 

Once we have a short chain of sufficient length (i.e. 
satisfying Eq. \A2\ we can attempt to join the chain to 
the large chain of length N. The expected number of 
qubits added to the long chain is thus mpgate — 1, at a 
cost of Rm + 1 attempted operations. The total cost is 
therefore 



C77 



Rr> 



1 



(A7) 



mp — 1 

attempted operations per qubit added to the long chain. 
Note that the dominant term in the denominator of this 
expression has an exponent that grows roughly loga- 
rithmically with Pgate"^- Thus, asymptotically, the cost 
grows slightly faster than polynomially with Pgate"^- 

The particular choice of m to be used depends on 
the value of Pgate- The reciprocal of Cm can be com- 
pared directly with the spillover rate in our numerical 
calculations, and therefore in Fig. [6] we plot the func- 
tions Ci7~^, Cg"^, Cs"^, C4~^ and Cs~^ for the regions 

1/17 < Pgate < 1/9, 1/9 < Pgate < 1/4, 1/4 < Pgate < 1/3 

and 1/3 < Pgate ^1/2, respectively. 

A similar calculation, which uses slightly more recy- 
cling and is therefore somewhat more efficient, was sub- 
sequently performed in Ref. [6] . Their calculations corre- 
spond to the non-deterministic CZ gate of Table O The 
approximate cost that they obtain for long linear chains 
can be determined from Eq. [2] of Ref. [g,] and is given by 



1 /oX log2(4/p+l) 



(A8) 



Again, C ^ can be compared directly with the spillover 
rate for the non-deterministic CZ gate, and is plotted in 
Fig.El 
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